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ABSTRACT 

(f) ' In this Paper we present a clustering analysis of QSOs as a function of luminosity 

over the redshift range z = 0.3 — 2.9. We use a sample of 10566 QSOs taken from 
the preliminary data release catalogue of the 2dF QSO Redshift Survey (2QZ). We 
analyse QSO clustering as a function of apparent magnitude. The strong luminosity 
evolution of QSOs means that this is approximately equivalent to analysing the data 
as a function of absolute magnitude relative to M* over the redshift range that the 
2QZ probes. Over the relatively narrow range in apparent magnitude of the 2QZ we 
find no significant (> 2er) variation in the strength of clustering, however, there is 
' marginal evidence for QSOs with brighter apparent magnitudes having a stronger 

, clustering amplitude. QSOs with 18.25 < bj < 19.80 show a correlation scale length 

so = 5.50 ± 0.79 h" 1 Mpc in an Einstein-de Sitter (EdS) universe and sq = 8.37 ± 
1.17 h^ 1 Mpc in a universe with £Iq = 0.3 and Ao = 0.7 (A), while the best fit values 
for the full magnitude interval (18.25 < 6j < 20.85) over the same spatial scales are 
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O ; s = 4.291^ h' 1 Mpc (EdS) and s = 6.35lg;g /i" 1 Mpc (A). We can therefore 

,H , determine that the bias of the brightest sub-sample is a factor 1.22 ± 0.15 (EdS) or 

1.24 ± 0.15 (A) larger than that of the full data set. An increase in clustering with 
luminosity, if confirmed, would be in qualitative agreement with models in which the 
luminosity of a QSO is correlated to the mass of the dark halo in which it resides, 
implying that the mass of the host plays at least some part in determining a QSO's 
r> \ formation and evolution. These models predict that the clustering in brighter QSO 

data sets, such as Sloan Digital Sky Survey QSO sample or the bright extension of 
the 2QZ should show a higher clustering amplitude than the 2QZ. 

Key words: galaxies: clustering - quasars: general - cosmology: observations - large- 
scale structure of Universe. 



1 INTRODUCTION 
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on hierarchical growth of structure predict that the level 
of clustering of a population is dependent on the mass of 
the dark halo in which the object resides. Thus clusters of 
galaxies should cluster much more strongly that galaxies do, 
which we observe to be the case. Analysis of galaxy surveys 
have also shown that more luminous (and therefore on av- 
erage more massive) galaxies have stronger clustering than 
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faint galaxies (Giavalisco & Dickinson 2001; Loveday et al 



The luminosity of QSOs appears to be (weakly) corre- 
lated with host galaxy luminosity, at least for low redshift 
QSOs (e.g. Schade, Boyle & Letawsky 2000). More luminous 
QSOs reside in brighter host galaxies. It is then natural to 
suppose that brighter QSOs should be found in more massive 
galaxies, assuming that luminosity and mass are correlated 
at least to some extent. Further circumstantial support for 
this argument is given by the correlation between the esti- 
mates of central blac k hole mass and the sp heroidal compo- 
nent of the galaxies ( Magorrian et al. 1998| ) . We would then 
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Figure 1. The redshift vs. absolute magnitude distribution of 
2QZ QSOs used in our analysis for an EdS cosmology. The solid 
line shows the value of M* as a function of redshift, as derived 
from the best fit LF model (see text), the model fit is calculated 
between z = 0.3 and 2.3, beyond z = 2.3 the model is an extrap- 
olation (dashed line). 



expect more luminous QSOs to cluster more strongly than 
faint QSOs. 

Until recently it was impossible to determine QSO clus- 
tering as a function of luminosity. The sparse nature of QSO 
surveys and the small numbers in homogeneous, complete 
samples meant t hat QSO clustering was only detectable at 



the ~ 4cr level ([iovino fc Shaver 1988 ; Andreani & Cris- 



tiani 19 ^3; |5hanks fe Boyle 1994 |Croom fe Shanks 199€| ; |La 
Franca et al. 1998|) . However, the 2dF QSO Redshift Sur- 



vey (2QZ) has allowed a dramatic improvement to be made 
in the measurement of QSO clustering. Using over 10000 
QSOs from the 2QZ, Croom et al. 2001a (henceforth Paper 
II) have made the first accurate (to ~ 10 per cent) mea- 
surement of QSO clustering. They find that QSO clustering 
averaged over the redshift range 0.3 < z < 2.9 is very simi- 
lar to that measured in redshift surveys of local [z ~ 0.05) 
galaxies. When fitting a standard power law of the form 

and 7 



(s/s )" 7 Croom et al. find s = 3.99Z 034 h' 1 Mpc 



l-58_o 09 f° r an Einstein-de Sitter Universe (hence- 
forth denoted by EdS). For a cosmology with Oo = 0.3 
and Ao = 0.7 (henceforth denoted by A) they find so = 
5.69±°:5o h~ x Mpc and 7 = 1.56±g;Jg- When investigated as 
a function of redshift the clustering of QSOs was found to 
be constant over the whole redshift range considered. 

The nature of any flux limited sample means that the 
most distant objects in the data set will on average have 
the highest intrinsic luminosity. There is then the reason- 
able concern that a comparison of QSO or galaxy clustering 
at different redshifts in any flux limited sample is actually 
measuring the properties of a different population of objects 
at each redshift. In the case of QSOs, their extreme evolu- 
tion in luminosity [oc (1 + z) 3 ] out to z ~ 2 (Boyle et al. 
2000; henceforth Paper I) means that at each redshift we 



are at least studying the same part of the QSO luminosity 
function (LF). 

In this paper we will make the first attempt to deter- 
mine the strength of QSO clustering as a function of lumi- 
nosity. In Section ^ we describe the data and analysis used. 
In Section^ we present our clustering results from the 2QZ, 
these are discussed in Section m 



2 DATA AND ANALYSIS 
2.1 The QSO sample 

For the analysis in this paper we have used the first pub- 
lic release catalogue of the 2QZ, the 10k catalogue (Croom 
et al. 2001b; henceforth Paper V). We only use QSOs with 
high quality (quality class 1; see Paper V) identifications of 
which there 10689 in the 10k catalogue (note that the num- 
ber of QSOs used in the analysis of Paper II was 10681, 8 
further QSOs were added to the 10k sample prior to its final 
publication). This 10k catalogue contains the most spectro- 
scopically complete fields observed prior to November 2000 
and is n ow publically available to the astronomical commu- 
nity at http://www.2dfquasar.org. The sample contains 



10566 QSOs in the redshift range 0.3 < z < 2.9 which will 
be included in our analysis below. 

The magnitude range of the 2QZ is 18.25 < fej < 20.85. 
This is a fairly narrow range, covering just over an order 
of magnitude in flux. However it does span the important 
regime in which the QSO luminosity function flattens to- 
wards fainter magnitudes. The distribution of QSO absolute 
magnitudes verses redshift is shown in Fig. [j]for the EdS cos- 
mology. We convert from apparent to absolute luminosity 
using the k-correction of Cristiani & Vio (1990) and assume 
a Hubble constant of Ho = 50 km s~ Mpc - (for luminosity 
determinations only) . We note that in our clustering analysis 
we include the factor h where Ho = 100ft -1 km s _1 Mpc -1 . 
The objects at high redshift are on average much more lumi- 
nous, as we would expect in a flux limited sample. Also, as 
the number of QSOs drops dramatically at brighter appar- 
ent luminosities there is only a small amount of overlap in 
intrinsic luminosity between QSO samples selected at widely 
differing redshifts. 

The QSO LF is well described by a double power law. 
The characteristic absolute magnitude of the break is de- 
noted by M* which is found to evolve strongly with red- 
shift. Using the power law evolution model of Paper I, in 
which M*(z) = M*(0) - 2.5(fci2 + k 2 z 2 ), Paper V de- 
rived the best fit parameters for the 2QZ 10k cat alogue 
and Large Bright Quasar Survey (Hewett et al. 1995) com- 



bined in an EdS universe. The parameters found were: $* = 
0.2 x 10~ 5 Mpc- 3 mag-\ a = 3.28, (3 = 1.08, Mb = -21.45, 
ki = 1.41 and k 2 = —0.29. M* as a function of redshift for 
this model is plotted in Fig. |l|. The fit to the LF is only 
carried out below z = 2.3, above that redshift the dashed 
line shows the extrapolation of the best fit model, which is 
approximately consistent with the density of high redshift 
( 3.6 < z < 5) QS Os found in the Sloan Digital Sky Survey 
( Fan et al. 2001 ). We have also fit the same model to the 
2QZ and LBQS assuming a A cosmology. In this case we 
find: $* = 0.46 x 10" 6 Mpc" 3 mag" 1 , a = 3.42, (3 = 1.35, 



-22.54, fei = 1.34 and k 2 



-0.26. This fit is in 
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Table 1. 2QZ clustering results as a function of limiting apparent magnitude. The fits assume a fixed power law slope equal to that of 
the best fit from the full data set. The listed \ 2 value is a reduced \ 2 f° r the best fit. 
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Figure 2. The QSO redshift distribution for our three apparent 
magnitude intervals, 18.25 < 6j < 19.80 (solid line), 19.80 < ftj < 
20.40 (dotted line) and 20.40 < 6j < 20.85 (dashed line). 



the range of acceptability specified in Paper I (Pks > 0.01) 
with a 2-D Kolmogorov-Smirnov (KS) probability of 0.02. It 
is worth noting that the parameters of these fits are strongly 
correlated, hence the use of any single parameter is to be dis- 
couraged. As this is the case we make comparisons between 
QSO samples of different apparent magnitudes, and do not 
split the QSOs in luminosity relative to M* . 

The division of QSOs on the basis of their apparent 
magnitude has several advantages over using absolute lu- 
minosities. First, apparent magnitude is a measured, not 
derived, quantity and is thus not dependent on the cosmo- 
logical model used. It should be easy to produce flux limited 
samples from any model of QSO formation, allowing direct 
comparison between data and models. Second, because M* 
is at approximately the same apparent magnitude at every 
redshift, an apparent magnitude cut selects the same part 
of the QSO LF at each redshift. 

Measuring the clustering of QSOs as a function of lumi- 
nosity is challenging. In particular, at luminosities brighter 
that M* the number density of QSOs falls off as a steep 
power law. Splitting a sample into N redshift bins of equal 



numbers of QSOs increases the error on pair counts by 
~ y/N. Dividing the sample into N luminosity bins increases 
errors by ~ NyN, as the surface density of sources is also 
reduced. This is a further reason why we limit our analysis in 
this paper to clustering as a function of apparent magnitude, 
as to carry out a meaningful sub-division on the basis of ab- 
solute magnitude requires a larger number of bins (due to 
the redshift- luminosity degeneracy inherent in the data set). 
A method to reduce the error in the measured clustering as 
a function of absolute luminosity is to cross-correlate QSOs 
of a particular luminosity with all other QSOs. This method 
will be investigated by Loaring et al. (in preparation). 

2.2 Correlation function estimates 

Our correlation function estimation is carried out as de- 
scribed in Paper II. For completeness we outline the pro- 
cedure here. The QSO correlation function, £q(s), where s 
is the redshift-space separation of two QSOs, is calculated 
using two representative cosmologies, the EdS and A models. 
The minimum variance estimator of Landy & Szalay (1993) 
is used to derive £q(s). In this Paper we divide the 2QZ 
10k sample into 3 apparent magnitude intervals with ap- 
proximately equal numbers of QSOs. The magnitude ranges, 
mean redshifts and other parameters for these bins are listed 
in Table [[]. We note that the mean redshifts of the differ- 
ent magnitude slices are very similar, with z — 1.398, 1.509 
and 1.552 for the bright, middle and faint slices respectively. 
These small differences cannot introduce a significant vari- 
ation in clustering, particularly as we find that QSO clus- 
tering does not evolve significantly with redshift (Paper II). 
The redshift distributions are shown in Fig. |^. 

We correct for the current incomplete observational cov- 
erage of the survey by using a random catalogue which ex- 
actly traces the distribution of observed QSOs on the sky, 
as in Paper II. The redshift distribution of these random 
points is taken from a spline fit to the QSO n(z) distribu- 
tion (each magnitude slice is fit separately). We also nor- 
malize the number of randoms within each UKST field in 
the survey to remove any possible systematic errors due to 
zero-point calibration errors. 

We calculate the errors on £q using the Poisson estimate 

of 



y/QQ(s 



(1) 



At small scales, <J50 h 1 Mpc, this estimate is accurate be- 
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Figure 3. The clustering of QSOs as a function of limiting apparent magnitude in an EdS universe, a) bright QSOs with 18.25 < 6j < 19.8, 
b) intermediate QSOs with 19.8 < fej < 20.4 and c) faint QSOs with 20.4 < 6j < 20.85. The solid line is the best fit power law in each 
case (assuming a slope of —1.58). The dashed line is shown for reference and is the best fit power law for the full sample from Paper II. 
The points at §(s) = 0.001 with no errors denote bins where no QSO-QSO pairs were found. These are properly taken into account in 
the fitting process. 



cause each QSO pair is independent (i.e. the QSOs are not 
generally part of another pair at scales smaller than this). 
On scales larger than ~ 50 h^ 1 Mpc the QSO pairs be- 
come more correlated and we use the approximation that 
A£q(s) = [1 + £q(s)]/ -J Nq, wher e No is the total num- 



ber of QSOs used in t he analysis (Shanks & Boyle 1994 



Croom fc Shanks 1996), for bins in which QQ(s) > Nq. 



Note that in this paper we are concerned with analysis on 
small scales (< 25 h Mpc), where the Poisson error es- 
timates are applicable. The yj Nq errors are only used for 
displaying our correlation functions. On very small scales the 
number of QSO-QSO pairs can be <10. In this case simple 
root-n errors (Eq. hi) do not give the correct upper and lower 



confidence limits for a Poisson distribution. We use the for- 
mulae of Gehrels (1986) to estimate the Poisson confidence 
intervals for one-sided 84% upper and lower bounds (cor- 
responding to lcr for Gaussian statistics). These errors are 
applied to our data for QQ(s) < 20. By this point root-n 
errors adequately describe the Poisson distribution. 



3 QSO CLUSTERING AS A FUNCTION OF 
LUMINOSITY 

In Fig. ^| we show the clustering of 2QZ QSOs for three 
different apparent magnitude intervals (see Table in an 
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Figure 4. The clustering of QSOs as a function of limiting apparent magnitude in a A universe, a) bright QSOs with 18.25 < bj < 19.8, 
b) intermediate QSOs with 19.8 < 6j < 20.4 and c) faint QSOs with 20.4 < 6j < 20.85. The solid line is the best fit power law in each 
case (assuming a slope of —1.56). The dashed line is shown for reference and is the best fit power law for the full sample from Paper II. 



EdS universe. The best fit to the total sample found in Paper 
II is shown as the dotted line for reference. On scales < 
20 tr 1 Mpc it appears that the brighter QSOs (Fig. |a) 
show marginally higher clustering that fainter QSOs (Figs. 
^3 and c). In order to determine the significance of this 
difference we make two measurements. We first determine 
the integrated correlation function, £q, within a radius s max , 
which is defined as 

CQ(Smax) = -3 / £ Q (x)x 2 dx. (2) 

Smax J 

We use s max = 20 Mpc as in Paper II. The results of 
this are shown in Table 0. While Cq(20) is greatest for the 
bright sample, it is not significantly larger than the values 
found for fainter QSOs. 



While £q is a robust statistic, it does not take into ac- 
count the the shape of £q. It is possible that a stronger 
constraint may be obtainable by fitting a functional form 
to the data. We fit a standard power law of the form 
£q(s) = (s/so) ' to our results. We make the assumption 
that the slope of the power law does not vary with lumi- 
nosity and fix the slope of the power law, to be the best fit 
from the full data set, that is, 7 = 1.58 (EdS) or 7 = 1.56 
(A). We use a maximum likelihood estimator based on the 
Poisson probability distribution function, so that 
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Figure 5. The QSO clustering scale length, so as a function 
of mean apparent magnitude in an EdS universe (filled circles) 
and a A universe (open circles). The best fit galaxy luminosity 
dependence models (see Section ^) for EdS (solid line) and A 
(dotted line) cosmologies are also shown. 



is the likelihood, where v is the observed number of QSO- 
QSO pairs, fj, is the expectation value for a given model and 
is the number of bins fitted. We fit the data with bins 
Alog(r) = 0.1, although we note that varying the bin size 
by a factor of two makes no noticeable difference to the resul- 
tant fit. In practice we minimize the function S = — 21n(L), 
and determine the errors from the distribution of AS, where 
AS is assumed to be distributed as x 2 ■ This procedure does 
not give us an absolute measurement of the goodness-of-fit 
for a particular model. We therefore also derive a value of 
X 2 for each model fit in order to confirm that it is a reason- 
able description of the data. We carry out the fit on scales 
0.7 < s < 20 ft" 1 Mpc, the scales chosen being the smallest 
scale containing a QSO-QSO pair and the scale at which £q 
appears to break from a power law. 

The best fit models are show in Fig. |§] (solid lines). 
As suggested by the £q measurements, the brightest sample 
does show stronger clustering than the other samples, how- 
ever the difference is only marginally significant. It should 
also be noted that the lowest clustering amplitude is found 
in the intermediate sample, although the fit to this data set 
is consistent with that of the faintest sample. 

In Fig. |^ we plot the clustering of QSOs in a A uni- 
verse. This shows the same feature as the EdS case, with 
the brightest QSOs being the most strongly clustered. The 
measured £q(20) values (Table ji]) also show this. The sig- 
nificance of the difference between the faintest and bright- 
est sub-samples is now larger, but still only at the ~ 1.7a 
level. We carry out the power law fit (this time between 
0.8 < s < 25 ft" 1 Mpc) as above (solid lines), again find- 
ing only a marginal difference between different magnitude 
intervals. 

In Fig. ^ we plot the fitted values of so as a function 
of mean apparent magnitude. In order to show a convincing 



trend as a function of magnitude we would like to span a 
much larger range in magnitude than is possible with the 
current sample. 



4 DISCUSSION 

Fitting the full data set over the same range of scales as 
above we find that s = 4.29± 'j^ ft" 1 Mpc (EdS) and 
so = 6.351±g|| ft" 1 Mpc (A). Thus the brightest third 
of QSOs has a clustering scale length which is a factor of 
1.28 ± 0.20 (EdS) or 1.32 ± 0.20 (A) larger than that of the 
full sample. As b/b* — (so/s^) 7 ' 2 this then implies that the 
ratio of the biases is 1.22T0.15 (EdS) or 1.24 ±0.15 (A). We 
therefore find only weak evidence that bright QSOs cluster 
more strongly. 

We should also consider comparing the above results to 
other measurements of luminosity dependent clustering. The 
most accurate measurement of this to date is by Norberg et 
al. (2001) using data from the 2dF Galaxy Redshift Survey. 
Norberg et al. find that the clustering of z < 0.3 galaxies is 
a weak function of luminosity fainter than L*, but a much 
stronger function of luminosity for galaxies brighter than 
L*. They find that the relation o g /&* = 0.85 + 0.15L g /L* 
well describes the luminosity dependence of galaxy bias, b s , 
relative to that found for Lg galaxies. If there was a sim- 
ple one-to-one relation between galaxy luminosity and QSO 
luminosity, we would be able to make a straightforward com- 
parison of QSO clustering to this relation. Unfortunately, we 
know that there is only a weak correla tion (with large dis - 
persion) between these two quantities (schade et al. 2000). 



We can, however, at least check for consistency between the 
galaxy and QSO clustering luminosity dependence. Under 
the assumption that galaxy luminosity is more closely re- 
lated to halo mass than QSO luminosity is, we would ex- 
pect the change in QSO clustering with luminosity to be 
no more rapid than the change in galaxy clustering with 
luminosity. Converting the above relation for galaxy clus- 
tering to a form directly applicable in our case we find 
(so/s,*) 7/2 = 0.85 + 0.15 x io-°- 4 < m - m *), where m* is the 
apparent magnitude corresponding to M* . The value of m* 
varies in the range 19.3 to 19.6 depending on the particu- 
lar best fit LF model used. We assume throughout that b 
is not a function of scale and allow So to be a single free 
parameter, while fixing m* [m* = 19.6 (EdS) or 19.3 (A)]. 
The best fit values are Sg = 4.78 ± 0.57 ft" 1 Mpc (EdS) or 
So = 7.20 ± 0.86 ft" Mpc (A) and lines corresponding to 
these are shown in Fig. H. We find that the above empirical 
model for galaxies is also consistent with the clustering of 
QSOs in the 2QZ. However, this is un-surprising given the 
large errors on the QSO clustering measurements. 

A number of authors have attempted to model the clus- 
tering of QSOs, based on the Press- Schechter (1974) fo rmal- 



ism for describing the evolution of dark matter haloes ( Mar- 



tini & Weinberg 2001; Haiman & Hui 2001). These make the 



assumption that the luminosity of QSOs is correlated with 
mass of dark matter halo in which they reside. Although 
there is reasonable evidence that the mass of compact ob- 
jects (presumably black holes) is correlated with the mass 
of the spheroidal component of the host galaxy (e.g. Magor- 
rian et al. 1998), a correlation between AGN luminosity and 
host luminosity and/or mass is less evident. The time-scale 
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of QSO activity is the main parameter that controls ampli- 
tude of clustering in these models, with time-scales of order 
10 6 yea rs being most consistent with t he results found in Pa- 
per II (Kauffmann & Haehnelt 2001) The main arguments 
are based on the number- density of galaxies that have gone 
through a QSO phase. Similar number- density arguments 
suggest that, as brighter QSOs are rarer, they should cluster 
more strongly. Our current analysis does not have sufficient 
signal-to-noise to clearly demonstrate this, and a more de- 
tailed investigation will have to await brighter QSO survey s 
such as the Sloan Digital Sky Survey ( ^Schneider et aL 200S ) 
and the bright extension to the 2QZ being carried out with 
FLAIR and the 6-degree Field system on the UK Schmidt 
Telescope. Further investigation is also possible using other 
techniques, such as the cross-correlation of QSOs of different 
luminosities. This will be investigated in the 2QZ by Loaring 
et al. (in preparation). 
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